Pranck-Condon blockade and giant Fano factors in transport through single molecules 
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We show that Franck- Condon physics leads to a significant current suppression at low bias volt- 
ages (termed Franck- Condon blockade) in transport through single molecules with strong coupling 
between electronic and vibrational degrees of freedom. Transport in this regime is characterized by 
remarkably large Fano factors (10^-10"^ for realistic parameters), which arise due to avalanche-like 
transport of electrons. Avalanches occur in a self-similar manner over a wide range of time scales, 
leading to power-law dependences of the current noise on frequency and vibrational relaxation rate. 

PACS numbers: 73.23.Hk, 72.70.-|-m, 73.63.-b, 81.07.Nb 



Introduction. — The Fano factor F, defined as the ratio 
S{lj — 0)/2e |/| of the zero-frequency current noise and 
the classical Schottky result for Poissonian shot noise , 
contains information about the charge of the current- 
carrying particles, and about quantum correlations be- 
tween them. Typically, the Pauli principle leads to sup- 
pressed Fano factors {F < 1) for fermionic carriers 
However, super-Poissonian noise has been discovered in 
systems with instabilities j^jj and, very recently, 
with dynamical spin blockade |5(. In the present pa- 
per, we report surprisingly large Fano factors in trans- 
port through single molecules, which originate from the 
coupling of electronic and vibrational degrees of freedom, 
but are unrelated to instabilities. 

Recently, single-molecule devices with vibrational de- 
grees of freedom (phonons) have received much attention 
both experimentally and theoretically 0,0, 0|- Here, 
we investigate the regime of strong electron-phonon cou- 
pling by current- voltage {IV) and noise calculations. Wc 
find that strong electron-phonon coupling leads to a sig- 
nificant current suppression at low bias voltages, which 
we term Franck- Condon blockade. The striking transport 
properties of this regime are further elucidated by com- 
puting the current fluctuations. We find that systems 
with weak vibrational relaxation can exhibit giant Fano 
factors of the order of 10^-10'^. This occurrence of giant 
Fano factors is due to avalanche- like transport of elec- 
trons interrupted by long times without charge transfer. 
These avalanches occur in a self-similar manner over a 
wide range of time scales, which is reflected in power- 
law behavior of the noise amplitude over many orders of 
magnitude in frequency and relaxation rate. 

Model. — Our results are obtained for a model of a 
molecule weakly coupled to metallic source and drain 
electrodes. As usual, we take the relaxation in the leads 
to be sufficiently fast so that their electron distributions 
have the form of Fermi functions. Transport through the 
molecule is assumed to be dominated by a single, spin- 
degenerate electronic level with energy e (measured with 
respect to the zero-bias Fermi energies of the leads) in the 
presence of one vibrational mode with frequency u)q. (A 
similar treatment for oscillations of the molecule's center 



of mass is straightforward [13.) This is described by the 



Hamiltonian H 
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leads and molecule. Here, Coulomb blockade is taken 
into account via the charging energy U . We focus on 
the regime of strong Coulomb blockade, [/ — > 00, appro- 
priate when eV,kBT <gC U. The operator d„ (dj.) an- 
nihilates (creates) an electron with spin projection a on 
the molecule, — '^o-'^o- denotes the corresponding 
occupation- number operator. Similarly, Capa {clp^) an- 
nihilates (creates) an electron in lead a {a — L,R) with 
momentum p and spin projection a. 

Vibrational excitations are annihilated (created) by b 
(6^). They couple to the electric charge on the molecule 
by the term ~ nd{b^ + b), which can be eliminated by 
a canonical transformation |^, |^|, leading to a renor- 
malization of the parameters s and [/, and of the lead- 
molecule coupling ta ta exp[— A(6^ From now on, 
we refer to the renormalized parameters as e and U. 

The coupling between molecule and leads is parameter- 
ized by the tunneling matrix elements 11 and Ir. Here, 
we consider the weak coupling regime so that the energy 
broadening 7 of molecular levels due to i?mix is small, 
i.e. 7 <C kBTjhioo, and a perturbative treatment for 
Hrn ix in the framework of rate equations is appropriate 
[l8j. For simplicity, we assume a symmetric device with 
t]^ = tfi = to and identical voltage drops of V/2 across 
each junction. (Generalizations are straightforward and 
do not lead to qualitatively new results.) 

The molecular occupation probabilities for n addi- 
tional electrons and q excited phonons on the molecule 
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are determined from the rate equations 

^± g' qi^q ^ q " q^q' 

n' ,g' 



(a) 
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Here, the coefficient W^^^j denotes the golden-rule rate 
for a transition \n,q) — > \n' ,q' ), consisting of the bare 
electronic tunneling rate 2TTpt^lh (with p the density 
of states of the leads), a lead Fermi function, possible 
spin factors, and a Franck-Condon (FC) matrix element, 
e.g. "l^. The last term in Eq. 10), containing 

the equilibrimn phonon distribution Pq'^^ describes re- 
laxation of the phonon excitations with relaxation time 
r. Within the rate-equations approach, the current can 
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be written as / = J2n , 
where the additional subindex L in the rates denotes 
transitions occurring between the left lead and the 
molecule only. 

Franck-Condon blockade. — The strong dependence of 
the IVs on the electron-phonon coupling strength A is 
shown in Fig. ^a) for A = 1 (intermediate coupling) 
and A = 4 (strong coupling), as obtained from the rate- 
equation approach. We first consider the case of e = 0, 
i.e. the molecular single-particle level and the lead Fermi 
energies are aligned at zero bias. Then, for A = 1, the 
current increases sharply due to resonant tunneling when 
switching on a small bias voltage, and it exhibits the 
characteristic FC steps. In contrast, for A = 4 the current 
is significantly suppressed at low bias voltages. 

The current suppression originates from the behav- 
ior of the FC matrix elements determining the rates of 
phononic transitions qi — > 92- Their dependence on the 
electron-phonon coupling strength can be understood in 
terms of the overlap of two displaced harmonic oscilla- 
tor wave functions. The parameter A determines the 
magnitude of this displacement in units of the oscilla- 
tor length. For weak coupling, A 1, transitions mainly 
occur along the diagonal qi — > gi . For intermediate cou- 
pling, A « 1, the distribution of transition rates becomes 
wider, and transitions slightly off-diagonal are favored. 
For strong electron-phonon coupling, A ^ 1, the distri- 
bution widens considerably and a gap of exponentially 
suppressed transitions between low-lying phonon states 
opens, see Fig.^b). It is interesting to note that these in- 
gredients are also crucial in the context of phonon broad- 
ening of resonant-tunneling line shapes 

For temperatures A;bT <C fioJo, a tunneling event can 
increase the number of excited phonons by at most Aq — 
int(0.5e \V\ /hujo) due to energy conservation. Thus, for 
strong electron-phonon coupling and at low bias voltages, 
the system is trapped in a region of exponentially small 
transition rates. For equilibrated phonons (t = 0) this 
suppression dominates the IV until the bias voltage is 
high enough {eV ~ X^Hujq) to escape from the blockade 
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FIG. 1: (Color online) (a) Current-voltage characteristics for 
intermediate (A = 1) and strong (A = 4) electron-phonon cou- 
phng for e = and ksT = O.OSftojo with (eq.) and without 
(un.) relaxation of phonons. We find that strong electron- 
phonon couphng leads to a significant current suppression at 
low bias voltages. This Franck-Condon blockade arises from 
the behavior of the Franck-Condon rates for phonon transi- 
tions gi §2 plotted in (b). The rates Fg^gj, shown for 
A = 1 (left) and A = 4 (right), are given in units of 2-Kpt'^/h. 
For strong electron-phonon coupling, transitions between low 
lying phonon states are exponentially suppressed. The corre- 
sponding current suppression cannot be lifted by a gate volt- 
age, which may serve as a fingerprint of FC blockade. This is 
depicted in the plot of dl /dV in the V-Vg plane for unequi- 
librated phonons with A = 4 (c). The case of intermediate 
coupling with A = 1 (d) is shown for comparison. 



regime by transitions from zero phonons to highly excited 
phonon states. For unequilibrated phonons (r = 00) the 
blockade is less rigorous, since a tunneling event can leave 
the molecule in an excited phonon state and subsequent 
tunneling events can increase this excitation even further. 

There exist two experimental fingerprints of the FC 
blockade, which can be readily used to distinguish it from 
other low-bias current suppressions, such as off-resonance 
tunneling for e 7^ and Coulomb blockade. In ungated 
devices, the succession of FC step heights yields a finger- 
print of the origin of the blockade regime. For symmetric 
devices with strong phonon relaxation, the presence of 
FC blockade is reflected in increasing step heights when 
leaving the voltage range of current suppression. For 
asymmetric devices or weak phonon relaxation, this cri- 
terion does not hold for each single step but still gives a 
valid tendency for the step-height succession. (We point 
out that systems with strong electron-phonon coupling 
and weak relaxation can even exhibit NDR phenomena 
as indicated by negative values for dl/dV in Fig. ^c) in 
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FIG. 2: Zero-frequency (excess) Fano factor as a function 
of bias voltage with A = 4 for e = and fcsT = Q.O^Tiujo. 
Results are shown for equilibrated phonons (r = 0) and un- 
equilibrated phonons (r = oo). We find that the combined 
effect of strong electron-phonon coupling and weak phonon 
relaxation leads to giant Fano factors in the Franck-Condon 
blockade regime. 



red color.) When working with gated devices, a clear dis- 
tinction can be achieved by measuring dl/dV as a func- 
tion of both bias and gate voltage. Here, the FC blockade 
can be identified by an extended blockaded region in the 
V-Vg plane, whereas an e 7^ suppression and Coulomb 
blockade can always be lifted by an appropriate adjust- 
ment of the gate voltage, see Figs.nic),(d). 

Giant Fano factors. — The remarkable transport prop- 
erties of the system in the FC blockade regime are elu- 
cidated by computing the current fluctuations. We use 
Monte Carlo (MC) methods to simulate the explicit tem- 
poral dynamics of tunneling events, and rate equations 
combined with a Langevin approach [l^ to calculate the 
current noise 



dte"^'{5I{t + t')5I{t'))t 



(3) 



Here, we restrict our discussion to symmetric junctions 
so that according to the Ramo-Shockley theorem the 
current is given by I{t) = {Iiit) + ^fi(^))/2. Typical re- 
sults for strong and weak phonon relaxation are depicted 
in Fig. 12 where the zero-frequency Fano factor Fex for 
the excess noise S — S\v=o is plotted as a function of bias 
voltage. For equilibrated phonons, Fex is suppressed be- 
low 1 and except for weak step-like structures, it remains 
essentially constant in close vicinity of F = 5/9 expected 
for the phononless system with one spin-degenerate level 
and strong Coulomb blockade 0|. 

In contrast, for unequilibrated phonons with A = 4 the 
Fano factor reaches values close to 200 as soon as the 
bias voltage is high enough for exciting the first phonon, 
and subsequently decreases stepwise at bias voltages cor- 
responding to the opening of further phonon channels. 
(We note that we found even larger Fano factors of the 
order of 10'^ for A = 4 and nonzero single-particle en- 
ergy e 7^ 0.) For bias voltages sufficient to lift the FC 
blockade, the Fano factor returns to values of the order 
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FIG. 3: Number of transmitted electrons A'^ and number of ex- 
cited phonons q as functions of time, obtained by MC simula- 
tion for T = 0.05fttJo, V = 3/itJo/e, e = 0, r = 00, and A = 4. 
The upper plot shows a magnification of the marked region in 
the lower plot. The results show long periods of time without 
current flow when the system assumes the phononic ground 
state, and avalanche-like current flow accompanied by phonon 
excitations on much shorter time scales. Self-similarity of the 
curves is found when magnifying the avalanche phases. 



of 1. 

The physics of the giant Fano factor originates from 
the combined effects of FC blockade and weak phonon 
relaxation, and can be understood by considering the 
time-dependent dynamics of the system. Results from 
a MC simulation of the tunneling events for ksT <^ Hcuo 
and strong electron-phonon coupling are shown in Fig.|2| 
(bottom), where the number of transmitted electrons N 
within the time interval [0, t] and the number of excited 
phonons q are plotted as functions of time t. We find 
that for large periods of time the system assumes the 
phononic ground state q = and the current vanishes 
(i.e. N is constant). These phases are interrupted by 
avalanche-like increases of N (accompanied by phonon 
excitations) on much shorter time scales. When mag- 
nifying the avalanche phases, we detect self-similarity of 
the curves N{t) and q{t), see Fig. 13 (top): Again, one ob- 
serves long times with constant N, now associated with 
(7 = 1, and short avalanche phases. 

The extended phases without charge transfer are due 
to the suppressed transition rates for low lying phonon 
states in the FC blockade regime. Typically, the system 
spends a long time in the phononic ground state before 
a transition takes place. Starting from q — and for 
voltages 2huJo < e\V\ < AHluq, where the Fano factor 
is maximal, energy conservation only allows the phonon 
transitions 0^0 and — > 1. Since the rates for phonon 
transitions — > 92 roughly grow exponentially with q2, 
it is favorable to increase the phonon number whenever 
an electronic transition occurs. Thus, as soon as the sys- 
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FIG. 4: Noise power in the Franck-Condon blockade regime 
as a function of (a) frequency, and (b) phonon relaxation rate 
for e = 0, A = 4, and ksT = O.OS^tJo- In (a) the total 
current noise is plotted for unequilibrated phonons. In (b) 
the zero-frequency noise is plotted for the bias voltage eV = 
ZhujQ. We find (approximate) power-law behavior S{uj) ~ 
w"" and 5'(0) ~ r^^" with a « 1/2, arising from the self- 
similar avalanches over a wide range of time scales. 

tern succeeds in leaving the phononic ground state and if 
phonon relaxation is weak, subsequent transitions tend to 
increase the phonon number even further, while the tran- 
sition rates increase roughly exponentially. Then, after 
having left the regime of suppressed transition rates, the 
system can exhibit an avalanche of tunneling events, un- 
til a transition to the phononic ground state induces an- 
other long period without charge transfer, and the cycle 
repeats. Effectively, this leads to electron bunching, and 
the Fano factor can be estimated by the typical number 
of electrons transmitted during an avalanche, see Fig. |21 

The remarkable self-similarity of avalanche phases due 
to the exponential stretching of transition rates is re- 
flected in the relaxation-rate dependence of the Fano fac- 
tor and the noise power spectrum, depicted in Fig. |3 
In both cases, we find approximate power-law behavior, 
S uj~°' and S ~ t^/°^ with a ~ 1/2, over many or- 
ders of magnitude in relaxation rate and frequency, re- 
spectively. (In some cases with A = 4 and e ^ 0, we 
found this power-law behavior to extend over 6 orders of 
magnitude in oj and 1/t.) The boundaries of the power- 
law scaling in the frequency domain are connected to the 
maximum available transition rate, and to the "smallest" 
transition rate describing the escape from the phononic 
ground state. Presumably, the weak oscillatory devia- 
tions from a pure power law arise from the discreteness 
of the variable q. 

Conclusions. — Vibrational degrees of freedom can play 
an important role in transport through single molecules, 
and lead to transport properties very different from those 
of conventional semiconductor or metal nanostructures. 
In particular, we have shown that FC physics, which is 
characteristic of molecules, leads to a low-bias current 



suppression (FC blockade) and giant Fano factors of the 
order of 10^-10'^ in single-molecule devices with strong 
electron-phonon coupling and weak phonon relaxation. 
The Fano factor enhancement has been explained by 
avalanche-like transport of electrons. The self-similar oc- 
currence of avalanches over a wide range of time scales is 
reflected in approximate power-law behavior of the noise 
as a function of frequency and relaxation rate. We fi- 
nally remark that an experimental verification of this new 
tranport regime depends on designing a single-molecule 
device with strong electron-phonon coupling and weak di- 
rect vibrational relaxation. Strong electron-phonon cou- 
pling is even realized in simple diatomic molecules such 
as F2 and Kr2 for which we estimate A ~ 4.4 and A ~ 5.4 
|l4l |. and very likely in a plethora of larger molecules. 
In addition, vibrational relaxation times as large as 10ns 
have been observed in recent experiments |l5j . For these 
relaxation times, the regime discussed in this paper is 
reached for currents large compared to lOpA [l^. Al- 
ternatively, the effects predicted here may be relevant 
in artificial nano-electromechanical devices. Indeed, the 
current suppression in the FC blockade may have recently 
been observed in such systems 
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